Spontaneous symmetry breaking and Lifshitz transition in bilayer graphene 
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We derive the renormalization group equations describing all the short-range interactions in bilayer 
graphene allowed by symmetry and the long range Coulomb interaction. For certain range of 
parameters, we predict the first order phase transition to the uniaxially deformed gapless state 
accompanied by the change of the topology of the electron spectrum. 
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Introduction- The Lifshitz transition (LiTr) [1] is the 
simplest topological effect in physics of metals. It con- 
sists of the change of connectivity of isoenergetic surfaces, 
either as a function of electron density or external param- 
eters, such as strain. As the change of the topology of 
the e.g. Fermi surface can not be continuous, all the ob- 
servables in the system should experience singularities at 
the LiTr also known as a half-integer-order phase transi- 
tion (PT) . Alternatively, the reconstruction of the Fermi 
surface may occur via an underlaying spontaneous sym- 
metry breaking PT. The observation of the LiTr in the 
bulk metals is an extremely challenging task: a varia- 
tion of the Fermi level in metals requires doping which 
introduces disorder and obscures the transition, whereas 
application of strain requires high pressure experiments. 




FIG. 1. Constant energy lines for the one particle spectrum 
in the BLG. LiTr as a function of density occurs when the 
Fermi level intersects the separatrix (bold line) for a) un- 
broken Csv symmetry of the graphene lattice; b) symmetry 
breaking Cav — >■ C2v. Circles mark Dirac points and ± indi- 
cates Berry phases, (pB = ±7r. Insets shows the tight-binding 
cartoons for the band structure (i) and the schematic evo- 
lutions of the Landau levels in the magnetic field (ii) (the 
number indicates the degeneracy per one spin and valley). 

Bilayer graphene (BLG) - a two-dimensional allotrope 
of carbon with a honeycomb lattice - is a potentially ideal 
system to study the LiTr [2] . A gapless low-energy elec- 



tronic structure of the conduction and valence bands near 
the Brillouin zone (Bz) corners in Bernal stacked BLG 
has a parabolic dispersion e± w ±p^/27n at intermedi- 
ate energies determined by the intra- and interlayer hops 
between closest neighbors. Remarkably, the electronic 
wavefunctions accumulate the ~ 2n Berry phase as 
the momentum going along the loop encompassing p = 0. 
This causes the double degeneracy (per one spin and one 
valley) of the zero-energy Landau level (LL) in the mag- 
netic field. Those features, however, are not protected by 
the crystal symmetry. The parabolic dispersion is trig- 
onally deformed at the lowest energies and e = state 
splits into four Dirac points: one in the Bz corner and 
three off-sets separated by momentum 2mw3 due to the 
next-neighbor interlayer hopping. Fig. la (this separation 
is about 0.1% of the size of Bz). The total is con- 
served so that the Dirac point in the Bz corner carries 
(j)B = — TT and each of the off-sets (j)B = tt. This doubles 
the degeneracy of zero-energy LL at weak magnetic field. 

Unlike in conventional metals, the LiTr in the BLG 
(from one to four Fermi lines) can be tuned by a small 
density variation controlled by a gate voltage. Suitable 
suspended BLG devices of sufficient quality for the LiTr 
studies have been fabricated [3, 4]. 

The central question of this Letter in the stability of 
the above Lifshitz transition against the effect of electron- 
electron interaction (EEI). There are two possibilities: 
(i) EEI does not break the symmetry leading to a quan- 
titative renormalization of the band structure affecting, 
e.g. the density tlutt = (2/7r^)(mD3/?i)^ corresponding 
to the LiTr; (ii) EEI does break the symmetry leading 
to a qualitative transformation of the spectrum - the 
number of the Dirac points is then determined by the re- 
duced symmetry (contradicting scenarios were suggested 
in Refs. [5-7] ). Using the renormalization group (RG) 
treatment of the EEI problem we found that: (i) n^iTr is 
not renormalized; (ii) the most likely spontaneous sym- 
metry breaking in BLG occurs by the generation of the 
asymmetric hopping in the effective Hamiltonian with 
the same symmetry as the effect of A — B sublattice dis- 
placement [7], see Fig. lb. For the gedanken experi- 
ment where vs is varied, the symmetry breaking occurs 
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via first order quantum PT after which the spectrum re- 
mains gapless but two Dirac points are annihilated and 
two other persist and carry (f>B = tJ"; as the result, the 
degeneracy of the zero-energy Landau level is half that 
of the for unbroken Cgy symmetry. The Berry phases 
control the degeneracy of the zero-energy LL which is 
visible via the Shubnikov-de Haas oscillations, and for 
some range of parameters the finite temperature PT is of 
the first order leading to bistabilities in transport. 

The low energy model for the bilayer graphene is for- 
mulated in terms of the states close to K and K' points 
of the Brillouin zone [2] . The Hamiltonian is 



H= cfr^l ho + hui + hc + hs 



(la) 



Hereinafter, the summation over repeated spin indices 
a = -ill 12 is implied. Four component fermionic field 

{KK') and the sublattice {AB) spaces [2] (sublattices A 
and B belong to the different layers) . All matrices acting 
in this four dimensional space are represented as direct 
product of the Pauli matrices f/^^,f/^^ , (1=0,1,2,3): 



where subscripts indicate the irreducible representations 
of the extended point group, see e.g. Sec. Ill of Ref. [9]. 
For example, Ex is two dimensional representation which 
does not change sign under C2 rotation and describes the 
symmetry breaking shown on Fig. lb, whereas B2 is the 
one dimensional representation describing the breaking 
of the interlayer symmetry, Cgt, — >■ C^^ . 

The RG study of the model (1) is based upon the analy- 
sis of diagrams shown in Fig. 2. The Coulomb interaction 
apparently is the most relevant operator (i.e. its pertur- 
bative treatment leads to the linear rather than the loga- 
rithmic divergence). The screening of this interaction, see 
Fig. 2c, makes it marginal; its value is ~ l/[iVn(g,a;)], 
(n is the polarization operator). The formal justifica- 
tion for the approximation Fig. 2c is the 1/iV expansion 
- which we believe is applicable for iV = 4 - and the long 
wavelength limit. Note, that enters together with the 
Coulomb interaction potential so that it drops out, see 
Fig. 2c, and does not contribute to the running of the 
coupling constants. Other constants gl are assumed to 
be small and treated in a first loop approximation [10]. 



= G(f-.k) = 



O =;io(t,fcj,); -O =L{k^,ky); 



and Tq ■ is the unit 2x2 matrix. 

The kinetic energy is given by (h = 1, k^. , 



(lb) 
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2 JV/g ky 



/(2m). (Ic) 



Together with Eq. (Ic), the trigonal warping term. 
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(Id) 



determines the spectrum in Fig. la. 
The long-range Coulomb interaction. 



e- I d\'^l,{v')^a,(Y') 
|r-r'| 



(le) 



is the strongest in the system. However, due to the 
screening it does not scale and therefore does not describe 
any symmetry breaking by itself. The latter is captured 
by the scaling of the marginal short-range interaction 




h,r = (2^/m) 9lMl y^^Ml^ 

ij=0 



(If) 



The couplings gf are not independent [8]. The Cgu sym- 
metry of the bilayer constrains 



FIG. 2. Derivation of the RG equations, a) Definitions of 
the elements; b) Polarization operator; c) Screening of the 
Coulomb interaction; d,e) Renormalization of one particle 
spectrum; f) gauge invariance of the scalar vertex; g) renor- 
malization of the short range interaction 
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Because the polarization operator does not have loga- 
rithmic divergences and all of the interactions are consid- 
ered in the first loop, the details of the cut-off scheme are 
not important. On each step we will restrict the inter- 
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nal momentum of the loop as £ — d£ < k'^/2m{£) < £. 
We, then, rescale '0 — ?> (1 + SZ/2)ip, {SZ is defined on 
Fig. 2d) to keep the term drip in the Matsubara equation 
intact. As a bonus, the scalar vertex is also not renormal- 
ized for the reason of gauge invariance, see Fig. 2f. Renor- 
malizations of m and V3, see Fig. 2e, are given by [11] 



dlnm/di ~ —dlnv^/di — —ai/N; ai 



.078, (2a) 



where i = \og{£o/£), where £q ~ 0.3 limits the appli- 
cability of the two-band model of bilayer graphene. 

The possible symmetry breakings are described by the 
scaling of the short range interaction terms [11]: 



dgf ^ a5{Ei)l _ aigl 
di iV2 N 



k.l.m.n—0 



gj = gi(l - S,oSjo) + 6,06,0^2/ {2N), as « .469 (2b) 

where a = — a^/lG, w .066, the symbol S(Ei)l 
is defined as S{Eiy-^l''^ = 1 and (5(£'i)^ — otherwise. 
The summation over repeated indices is not implied in 
Eq. (2b). The constants in Eq. (2b) are given by 
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Note that, Eqs. (2b) and (2c) respect symmetry (Ig). 

Equations (2) are the main technical result of this pa- 
per. They describe the evolution of all the band struc- 
ture parameters and all short-range interactions terms 
allowed by symmetry in the leading logarithmic approxi- 
mation. To compare with the existing literature: the RG 
treatment of Ref. [7] considers only two possible terms 
{gE^^ and gsi ), treats the Coulomb interaction as short- 
range and neglects the warping in the spectrum; mean- 
field treatment of Refs. [5, 6] corresponds to hardly jus- 
tifiable taking into account only one [5] [(i) of Fig. 2g[ 
or two [6[ diagrams [(i,ii) of Fig. 2g[ with the subsequent 
projection on the B2 representation. 

RG flow and non-broken symmetry - The density of 
electrons (or holes) at which the topology of the Fermi 
surface changes is found from Eqs. (Ic) and (Id) as 
riLiTr = (2/7r2) {mvz/fif ~ 2 x 10^°cm~'^ (estimated 
with m = 0.035 and = W3 ~ lO^cm/s ). According to 
Eqs. (2a), the Coulomb part of the EEI does not renor- 
malize n^iTr but affects the energy of the saddle points 
in the single-particle spectrum £LiTr = mv'^/2. The 
bare value of this energy can be estimated using the bi- 
layer parameters m, V3 quoted above as £LiTr — 1 meV. 



The renormalized value is £LiTr = £LiTr {£o/£LiTr) " — 
£LiTr {£LiTr/£o) , such chaugc is uot obscrvablc. 

RG flow and symmetry breaking - The divergence of 
a coupling constant g. during the renormalization sig- 
nals the symmetry breaking with the order parame- 
ter from the corresponding irreducible representation (a 
more complete classification, involving the the magnetic 
and gauge symmetries will be reported elsewhere [12]). 

Let us assume with the short-range interactions on the 
energy scale Eq are negligible, 5^ = 0. The constant 
term in Eq. (2b) means that this point is not fixed and 
couplings gEi and gsi (see Eq. (Ig)) will fiow away from 
this point. Ignoring gsi, we obtain an equation for gE^'- 

^-iv(FT2)-^<« + ^'<'"'.-«''^ (3) 
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Note that Ci > 0, for > 0, and no fixed point exists; 
though aa appears small, its neglecting would lead to a 
nontrivial fixed point g^^ ~ Solution of Eq. (3) is 



gEi (i) = C2 



Cl 



2N{N + 2) 



■ cot 



(4) 



where £0 is found from gEi{£ = 0) = 0: for A^ = 4 
£0 ^ 7.1. Inclusion of gs^ shifts the pole slightly so that 
(gE,; gsJoc(4-^)-H-l-67;0.85). 

This divergence implies a symmetry breaking at £ei — 
Eoe^'^^^ ~ 0.3meV [13[. It is important to notice that £ei 
and £LiTr turn out to be of the same order and, therefore, 
have to be considered together. A more accurate theoret- 
ical comparison of those two energy scales requires more 
detailed knowledge about the microscopic values of the 
initial interaction constants which is not available at this 
time. Therefore, we will discuss the possible PTs for an 
arbitrary value of T = £eiI £LiTr- 

Analysis 0/ the phase transition- If T < 1. the di- 
vergence of gEi is terminated and the symmetry is not 
broken. The possible divergence of the coupling constant 
gEi at T > 1 indicates the symmetry breaking and ap- 
pearance of the anomalous averages comprising the irre- 
ducible representation Ei of the group Cq^ , 



u, = {2TT/m){^lMi^^); j = l,2. 



(5) 



For studying the PT we have to consider the Landau free 
energy density. It must be of the form 



/ = nLiTr£LiTr^T 



'^LiTr 



£^ 

'-'LiTr 



(6) 



for the symmetry and dimensionality reasons. At T < Tc 
(here Tc ~ 1), function T{x,0) has a local minimum at 
2; = which, at T > Tc, turns to a maximum. The pres- 
ence of the cubic invariant prescribed by Cgt, symmetry 
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FIG. 3. a,b)The curvature of the mean-field energy along 
the steepest decent direction. The square root singularity is 
caused by the collision of the Dirac points shown on the inset. 
At T < Ti < Tc extra local minimum and maximum are 
formed (a), such as lu™"^! > EuTr (b), i.e. only two Dirac 
points remain. The total dashed area equals to zero. c,d) 
The schematic phase diagrams for the finite temperature (c) 
and for quantum (d) (controlled by the gate voltage Vg) PTs. 
The insets show predicted hysteretic (or slow noise) behavior 
(dashed areas) of the conductivity for the corresponding paths 
on the phase diagram. 



signals that the zero temperature PT, under varying T, 
can be only of the first order and occurs at T < Tc. 

Now, we argue that the value of the order parameter 
in the ordered phase is such that the electron spectrum 
has two Dirac points, as in Fig. lb. In the mean-field 
approximation, the one particle Hamiltonian reads 



(7) 



At u\ = ~£LiTr,U2 = two Dlrac points collide and 
disappear and the band structure of Fig. lb is formed. 
At ui ~ 3£LiTnU2 = three Dirac points collide, and, 
once again, the spectrum with two Dirac points is formed. 
The mean-field energy density is given by 



/m_f(ui, M2) = (Tc - T) 
d^k 



mUiUi 
27r 



(8) 



e(k) - e(k)|„i_2=o 



UiUj 9^e(k) 



2 duidui 



Ul,2=0 



where e is the negative eigenvalue of H, see Eq. (7), and 
the summation over the repeated indices i,j = 1,2 is 
implied. The curvature of the energy density, found from 
Eq. (8), see Fig. 3, indicates that indeed ui formed during 
the PT transforms spectrum of Fig. la to that of Fig. lb. 

The corrections to the mean-field can not remove the 
singularity for the colliding Dirac points, as the Hamilto- 
nian (7) at low energies is protected by symmetry. There- 
fore, we believe, that our conclusion about the number 
of Dirac points in ordered and disordered phases is more 
general than the mean-field derivation. 



The conclusion about the first-order quantum phase 
transition have the important consequences for the finite 
temperature phase diagram, see Fig. 3c. At low tempera- 
tures T <C ShiTr, the transition remains of the first order 
up to some tricritical temperature Tt, and at Tt the tran- 
sition is continuous and belongs to 3 states Potts model 
universality class [7]. The quantum phase transition can 
be studied as the function of density n, and the phase 
diagram is on Fig. 3d. The number of Dirac poiirts in 
the ordered phase remains the same. 

In conclusion^ we investigated the interplay of the 
trigonal spectrum of the bare spectrum of the bilayer 
graphene with the electron-electron interaction. The de- 
rived RG equations allowed us to reveal the rich phase 
diagram [12] determined by the few (currently unknown) 
microscopic inputs. For a reasonably wide range of the 
initial conditions [13] we found Cgt, Civ symmetry 
breaking and connected it with the change of the topol- 
ogy of the single particle spectrum. We predicted the 
quantum phase transition of the first order as a func- 
tion of the electron density. Such a transition should be 
most readily observed in the hysteretic dependence of the 
conductivity on the gate controlled carrier density in the 
vicinity of the neutrality point. 
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